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I. INTRODUCTION 

From theoretical and experimental observations it is known that reaction-diffusion systems play an important role 
in the study of generic spatiotemporal behavior of nonequilibrium systems. Usually such models admit main contri- 
butions related to both local dynamics (chemical reactions type of birth- and-death processes) and mass transport. 
Novel experimental methods, such as field ion microscopy, scanning tunneling microscopy can be used to monitor 
chemical reactions on the metal surfaces with atomic resolution. 

In adsorption-desorption processes when material can be deposited from the gaseous phase such experimental 
methods allow one to investigate formation of clusters or islands of adsorbed molecules/atoms [1 J. Such islands 
can have linear size of nanometer range [2 j. In Refs.[3H7] it was experimentally shown that nano-patterns on solid 
surface and nano-islands in adsorbed mono-atomic layers can be organized. It was found that patterns on scales 
shorter than the diffusion length emerge from the interplay of reactions and lateral interactions between adsorbed 
particles. The adsorbate presence can modify the local crystallographic structures of the substrate's surface layer 
producing long range interactions between adsorbed atoms and their clusters (see for example, Refs.[8HT0]). It was 
observed experimentally that nanometer-sized vacancy islands can be organized in a perfect triangular lattice when 
a single monolayer of Ag was exposed on Ru(0001) surface at room temperature [7 j. Nanometer elongated islands 
was observed experimentally in Si/Si(100) [IT] , the same was found at deposition of Ge on Si [12], metallic elongated 
islands were observed at deposition of Cu on Pd(llO) [13 . It was shown that elongated adsorbate clusters are governed 
by formation of dimers and their reconstructions [14] representing nonequilibrium chemical reactions. 

It is well known that short-range transient patterns can be observed at initial stages of phase separation processes 
[T5] . These transient patterns can be stabilized by nonequilibrium chemical reactions leading to emergence of sta- 
tionary patterns [16]. In models of systems with adsorption and desorption processes additional chemical reactions 
should be introduced to freeze the patterns. It was shown previously, that adsorption and thermal desorption are 
"equilibrium reactions" which can not induce the formation of kinetic spatially modulated stationary phases [T71 [18] . 
A problem of formation of stationary microstructures in such systems with irreversible nonequilibrium chemical reac- 
tions was considered in Refs.pjJJ [20] . Properties of pattern formation in systems of adsorption-desorption type with 
dissipative dynamics were studied previously p~8j [2T] . Pattern formation in pure dissipative stochastic systems with 
multiplicative noise obeying fluctuation-dissipation relation was discussed in Refs. [22-25 . It was shown that such 
multiplicative noise can sustain stationary patterns of nano-size range in pure dissipative systems. 
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As far as the reactions are defined through chemical kinetics, an evolution of the field variable, let say coverage for 
adsorption/desorption systems, is governed by the reaction-diffusion equation. A typical deterministic equation is of 
the form 

ft* = /(*)- V ■ J, (1) 

where x = x(r,t) is the local coverage at surface defined as the quotient between the number of adsorbed particles 
in a cell of the surface and the fixed number of available sites in each cell, x < 1. The term f(x) stands for local 
dynamics and describes birth- and-death or adsorption-desorption processes; the flux J represents the mass transport. 

Most of the theoretical studies deal with the standard Fick law J = -DVx, where D is the diffusion constant. It 
is known that at / = the ordinary diffusion equation d t x = V • DVx has the unrealistic feature of infinitely fast 
(infinite) propagation. In such a case a fundamental solution 

x(r,t) = (47rL>fT 1/2 exp(-r 2 /Dt) 

means that for any small t at any large r the quantity x will be nonzero, though exponentially small. It leads to 
unphysical effect that particles move with infinite speed (more than sound speed in solids). The reason for this is 
lack of correlations of particle motion. To avoid such pathology the diffusion flux can be generalized by taking into 
account memory effects [26] 



M(t,t')DVx(r,t')dt' (2) 

described by the memory kernel M(£, t') = rj 1 exp(— \t— t'\/rj). The delay time tj is related to correlated (persistent) 
random walk. In the limit f(x) = one gets the finite propagation speed y ' D jrj. At tj — >• the asymptotic 
M(t,t f ) = S(t — t') leads to the Fick law J = —DVx with infinite propagation. As far as real systems (molecules, 
atoms) ha ve fin ite propagation speed one should use Eq.Q or an equivalent equation: rjd t J = —J — DVx [27] . 
Equations (TJ2) can be combined into the one reaction-Cattaneo equation of the form 



Tjd 2 tt x + (1 - Tjf'(x))d t x = f(x) + V • DVx, (3) 

where prime denotes derivative with respect to the argument. In the absence of reaction term (/ = 0) the corresponding 
telegraph equation has a solution of the form 



x(r,t) 



^ eXP 




where N = y/ ADtj , x — (D/rj)t 2 — r 2 , Jo(*) is the modified Bessel function. As was pointed out in Ref.[27] this 
equation has some restrictions: (i) it typically does not preserve positivity of the solution x(r,t); (ii) the damping 
coefficient must be positive, i.e., f'(x) < rj 1 . 

Therefore, one can use more realistic model given by Eq.([3|, where particles have finite speed at smaller time scales 
and approach diffusion motion on larger time scales. It follows that depending on the form of the reaction term f(x) 
Eq.Q admits oscillatory solutions [27 . An application of such formalism for phase separation processes study with 
finite rj allows one to describe pattern selection processes at early stages of decomposition in binary systems (see 
Refs.[28, 29 ) and oscillatory formation of the ordered phase in crystalline systems [30]. Oscillatory solutions in class 
of reaction-Cattaneo systems with fluctuating quantity rj were considered in Ref . [3T] . 

In this study we are aimed to describe dynamics of pattern formation and selection processes in a class of reaction- 
Cattaneo systems given by Eq.Q. Novelty of our approach is in studying oscillatory dynamics of pattern formation 
in such class of models with chemical reactions governed by adsorption/desorption processes. Following the formalism 
proposed in Refs.[20, 22j [23] we compare a behaviour of ordinary reaction-diffusion dissipative system and reaction- 
Cattenao system. It will be shown below that for the last class of models pattern selection processes are realized, 
stable patterns possible only if nonequilibrium chemical reactions are introduced. Studying behavior of islands size 
as clusters of dense phase we will show that averaged island size behaves itself in oscillatory manner. Considering 
formation of islands of adsorbed particles we discuss properties of the island size distribution during the system 
evolution. 

The paper is organized as follows. In Section II we propose the stochastic reaction- Cattenao model. The linear 
stability analysis and properties of pattern selection processes are presented in Section III. We discuss results of 
numerical simulation in Section IV. Finally, in Section V, we draw conclusions from our study. 
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II. MODEL 



Let us consider a model where only one class of particles is possible. Following Refs. [18-23, 32 one assumes that 
the particles can be adsrobed, desorbed, can diffuse and interact among themselves. Therefore, we introduce the 
scalar field describing dynamics of the local coverage at surface x(r,t) G [0,1]. The reaction term incorporating 
adsorption and desorption terms together with nonequilibrium chemical reactions is as follows: fo(x) = k a p(l — x) — 
kdxexp(U(r)/T) — k r x n . Here k a and kd are adsorption and desorption rates, respectively; p is the partial pressure of 
the gaseous phase; U (r) is the interaction potential. The last term corresponds to nonequilibrium chemical reaction 
of the order n > 1 with the rate constant k r . Usually it estimates islands/dimers formation or associative desorption 
[14]. In further consideration we put n = 2. 

The total flux J is a sum of both ordinary diffusion flux (— DqVx) and flow of adsorbate (—{Dq/T)x(1 — x)VU). 
Here the multiplier x(l — x) denotes that the flux is only possible to the (1 — x) free sites. Hence, the total flux is 

J = -DqVx - y-x(1 - x)VU\ (4) 

in an equivalent form one has 

* = -°oM( X ){^ + ±Vu], (5) 

where the Cahn mobility M(x) = x(l — x) is introduced. 

Formally the right hand side of Eq.([5| can be rewritten as follows: 

J = -A)M(x)V^, (6) 
where the total mesoscopic free energy functional is 

J 7 — «7"b + Pint- (?) 

The non-interacting part takes the form 

J*o = J dr [x(r) ln(x(r)) + (1 - x(r)) ln(l - x(r))] , (8) 

whereas Ti n t is governed by the interaction potential U which we assume in the standard form [18] 

U(r) = - J dr'u(r - r')x(r'), (9) 

where —u(r) is the binary attraction potential for two adsorbate particles separated by the distance r, it is of 
symmetrical form , i.e. f drr 2n+1 i^(r) = 0, n = 1, . . . , oo. 

Therefore, one can rewrite the total mesoscopic free energy as follows: 

T = J dr [x \n(x) + (1 - x) ln(l - x)} - J J drdr / x(r)^(r - r r )x(r f ). (10) 

Following Ref.[20j as a simple approximation for the interaction potential, we choose the Gaussian profile 

u(r) = exp(-r 2 /4r'), (11) 

where e is the interaction strength, tq is the interaction radius. Assuming that x does not vary significantly within 
the interaction radius, one can estimate 



J dr^(r - r>(r') ~ J &y'u(y - v') ^ ^— ^V n x(r). (12) 



Substituting Eq.|TT| into Eq.(12) up to terms of the 4-th order one gets J u(r)x(r)dr = 2ex, ^ J u{r)r 2 V 2 x(r)dr 
2er 2 ) V 2 x, ^ J ^(r)r^V 4 x(r)dr = erffl^x. 
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Therefore, using notation e = e/T one has 

^ J dru(r - r')x(r) ~ ex(r) + e(l + r 2 , V 2 ) 2 x(r). 
The total free energy functional takes the form [40 



ar + x\nx + (1 - x) ln(l - x) - -x(l + roV^a; 



Therefore, the total flux can be written as 

J = -D M(x)V 



loc 



s(l + rlV 2 ) 2 x 



Sx 



(13) 



(14) 



(15) 



here T\ oc is defined through the local part of the free energy density (first three terms in Eq.(14)). 

Next, it will be more convenient to measure time in units introduce the diffusion length L d = ^Do/k d and 
dimensionless adsorption rates a = k a p/k dl /3 = k r /k d . Therefore, the reaction term takes the form f(x) — a(l — x) — 
xe~ 2ex — fix 2 and the system is described by two length scales[4l], where ro <C L d . As far as real systems (molecules, 
atoms) have finite propagation speed one should take into account memory (correlation) effects, assuming 



3 = -L d J dt'M(t,t';rj)V 



5T 



Sx(r, V) ' 



(16) 



where A4(t, t') is the memory kernel. Taking it in the exponential decaying form A4(t, t') = Tj 1 M(x(r, t')) exp(— \t — 
t'\/rj), where tj is the flux relaxation time instead of one equation for the coverage we get a system of two equations: 



d t x = f(x) - L d V • J; 
rd t J = -J - L d M(x)V — , r = rjk d . 

At r j -)> the asymptotic M(t,t') = S(t - t') leads to the Fick law J = 
propagation. 

The equivalent equation for the coverage takes the form 



(17) 



-L,iM(x)V5J r t t/5x with an infinite 



rdt t x + j(x)d t x = <f(x,V), -y(x) = 1 -r/'(x); 

<p(x; V) = f{x) + L 2 d V ■ [Vx - eM(ar)(Va; + VC SH x)}, C SH = (1 + r 2 V 2 ) 2 . 



(18) 



Let us consider stationary homogeneous system states. In the deterministic limit in the absence of reaction, the 
system undergoes first order phase transition, where its stationary uniform states are given by a(l — x) = xexp(— 2ex). 
The critical point of this equilibrium phase transition x c = 1/2 is located at e c = 2 and a c = exp(— 2). A coexistence 
line of diluted (x < 1/2) and dense (x > 1/2) phases is given by the relation a = exp(— e). Introduction of the 
chemical reactions governed by the rate (3 shifts the whole phase diagram and shrinks the domain where the system 
is bistable. Here possible values for x related to the uniform states decreases whereas critical values for e become 
larger; if f3 grows then critical values for a decrease. The corresponding phase diagram is shown in Figjl] Outside the 
cusp the system is monostable: at small a (before the cusp) the system is in low density state xld, whereas at large 
a the high density state xhd is realized. In the cusp the system is bistable. Critical points and the related coverage 
are shown in Fig{T]3. Here starting from the fixed value for (3 and moving up to an intersection with the curve £*(/?), 
one gets the critical value £* as the corresponding ordinate in the left axis. Than, we move to the right hand side to 
an intersection with the curve £*(a*); the corresponding abscissa in the top correspond to critical value a*. Moving 
from this point to an intersection the curve x*(a*), one gets the critical value x* as the corresponding ordinate in 
the right axis. Acting in such a manner, we obtain all needed critical values for the system parameters. Analytical 
relation between a*, e*, x* and f3 are as follows: £* = (1 + V$ ~ 2x*)/(2x*(2 - #*)), a* = 2e* 2 x* 3 exp(-2e*x*), 
(3 = 2e*(l - £*£*) exp(-2e*x*). 



III. LINEAR STABILITY ANALYSIS 



It is known that systems with memory effects admit pattern selection processes at fixed set of the system parameters 
[28H3Q] . These processes can be observed at early stages of the system evolution where linear effects are essential. 
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Therefore, pattern selection can be studied considering stability of statistical moments, reduced to the averaged filed 
and/or structure function as the Fourier transform of a two-point correlation function for the coverage. As far as 
fourth order contribution in the interaction potential u(r) is not essential at small ro < L^, next we consider a case 
where J u(r)x(r)dr ~ 2e(l + Tq\7 2 )x, neglecting TqV 4 x. 



Averaging Eq.(18) over initial conditions and taking (x) —xqoc e 1 ^ 1 kr ^ one gets the dispersion relation of the form 

1/2 



^W=f = 



L 2 P(1 - 2£M(x )(l - rgk 2 )) - /' (x ) 7 2 K) 



4r 2 



(19) 



One can see that uj(k) can have real and imaginary parts, i.e. uo(k) = ^iuj(k) ± i9u;(fc). The component SRu;(fc) is 
responsible for oscillatory solutions, whereas 9u;(fc) describes stability of the solution (5x\z(uj)}. Analysis of both 
$toj(k) and 9a; (k) allows us to set a threshold for a wave-number where oscillatory solutions are possible. Moreover, 
it gives a wave-number for the first unstable solution. From the obtained dispersion relations it follows that at k = ko 
satisfying equation 



^ 2 (l-2eM(xo)(l-r 2 ^ 2 )) = ^ 



/'(so) 



7 Oo) 



(20) 



4r 

two branches of the dispersion relation degenerate. The unstable mode appears at k = k c obtained from the equation 

L\k\(\ - 2eM(x )(l - rlkl)) = f(x ). (21) 

From the dispersion relation uo(k) one can find the most unstable mode k m as a solution of the equation d9o;/d/c = 0. 
It coincides with first unstable mode when only one nonzeroth solution of the equation 9a; (k) = emerges. 

Using obtained relations one can calculate a diagram indicating spatial stability of all homogeneous states to 
inhomogeneouse perturbations. The corresponding diagram is shown in Fig(2] Here domain of unstable modes with 
respect to inhomogeneous perturbations is limited by solid and dashed thick curves. The solid curve relates to high 
density phase, whereas dashed line corresponds to low density phase; dotted line addresses to unstable homogeneous 
stationary state. When we increase the adsorption rate a from zeroth value the first unstable mode emerges at large 
k and is possible only for high density phase. There is small domain for a where spatial instability of the low density 
phase is possible (see magnified insertion of a(k) at small a and k). It is should be noted that wave numbers related 
to these unstable modes in both low- and high density phases are observed in fixed interval k G [fc c i,fc C 2]. The thin 
solid line in Fig j5] denotes critical values k Q where oscillatory solutions (x(k,i)) are possible. The domain of unstable 
modes with respect to inhomogeneous perturbations is limited by large values for a. It means that instability of high 
density phase is possible only in fixed interval for adsorption rate values. 

According to obtained dependencies a(k) we plot in Fig J3] critical values for a and e related to formation of spatially 
modulated phases, where solid lines denote binodals, dashed lines bound spatially modulated phases. In the case of 
/3 = one has only binodals bounding domains of one uniform phase and domain of two uniform phase existence. 
Here no spatially modulated phases are possible. When we put j3 ^ and increase the adsorption rate at small a 
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FIG. 3: Critical values for e and a related to pattern formation: a) f3 = 0.05; b) /3 — 0.1 



only uniform low density state (uLD) is possible; at large a one has only uniform high density state (uHD). In the 
cusp depending on the values for a one has: uniform low density state and modulated high density phase (uLD& 
mHD) at small a; modulated low density and high density phases (mLD& mHD) at elevated a. If we further increase 
the adsorption rate only modulated high density phase (mHD) is possible. Controlling the rate f3 one can govern size 
of the domain where both modulated low density and high density phases emerge; when f3 grows the domain of the 
modulated high density phase increases and mHD phase can emerge before the first binodal (at smaller a = a^i). 

Considering properties of pattern selection we need to find dynamical equation for the structure function 5(k, t) as 
the Fourier transform of the two-point correlation function (5x(r, t)5x(r' , t)) and study its behavior at small times. 
To that end we obtain the linearized evolution equation for the Fourier components Sx^(t) and Sx-\^(t) and compute 
S(k, t) = (Sx\^(t)Sx-\^(t)). The corresponding dynamical equation takes the form 

rd 2 t S(k,t) +1 (x )d t S(K t) = 2 {/'(x ) - L 2 d k\l - 2eM(x )(l - r 2 k 2 ))} 5(k,t). (22) 

Analytical solution can be found assuming S(k, t) — So ex exp(— izu(k)t), where 

2 (L 2 k 2 (l - 2eM(x )(l - rgfc 2 )) - f(x )) ^(x )] 1/2 



As in the previous case ^svj is responsible for stability of the system, whereas IStw relates to pattern selection processes. 
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x=0, e=4,a=0.1 x=0.5, 8=4, a=0.06 x=0.5, e=4, a=0.1 x=0.5, 8=3, a=0.1 x=0.5, 8=3, a=0.25 




FIG. 4: Snapshots of the system evolution at /3 = 0.1. All snapshots are taken at t = 0, t — 20, £ = 60, t — 160, t — 300. 
Here concentration of adatoms is shown with the help of gray scale: white domains correspond to adatoms present, dark ones 
indicate regions without adatoms. 

IV. NUMERICAL SIMULATIONS 

A morphology of emergent patterns and dynamics of pattern formation we study by numerical simulations in two 
dimensional system with 256 x 256 sites and periodic boundary conditions. In our simulations we take time step 
At = 2.5 x 10 -4 and consider the case when — 40ro- The total size of the system is L = 12. SL^. As initial 
conditions we take: (x(r,0)) =0, ((Sx(r,0)) 2 ) =0.1 where x(r,t) G [0,1]. 

Typical evolution of the system with different values for r, e and a at /? = 0.1 is shown in Figj4j In the case of pure 
dissipative system (r = 0) at a = 0.1 the adsorbate is organized into separated islands with small difference in the 
linear sizes of islands. When we put r = 0.5 (see 3-rd column) at the same other system parameters we get pattern 
where all islands at large time interval are combined into one percolating cluster. Fixing r = 0.5 and taking small a 
one gets set of adsorbate islands with large difference in their linear sizes (see 2-nd column). With further increase in 
a islands of vacancies are formed (see 4-th column). Such islands has the same structure as islands of adsorbate at 
small adsorption rate. 

In Figure [5] we present dynamics of the averaged concentration field (x) and its dispersion ((Sx) 2 ). It is known that 
the growth of the quantity ((Sx) 2 ) means ordering of the system (an increase in fluctuations of the concentration). 
Following the obtained dependencies one can find that if nonequilibrium reactions are absent (/? = 0) the system 
passes toward an equilibrium thermodynamic state where no stationary patterns can be realized. Here all possible 
patterns appeared during the system evolution are transient and at final stages the system is totally homogeneous. 
According to the behavior of both (x) and ((Sx) 2 ) one can say that the average goes toward uniform stationary value 
xo — xhd, whereas the dispersion increases at early stages (formation of transient patterns) and after it decreases 
toward zeroth values (no dispersion in adatoms concentration field is realized). Therefore, the systems moves into 
homogeneous state. The situation is crucially changed when the nonequilibrium chemical reactions are introduced. 
Here such reactions freeze patterns formed at coarsening stage, and patterns develop very slowly. As a results of 
competing between chemical reactions and potential interactions between atoms such patterns are stable in time and 
can be considered as stationary ones (see Fig 13. Indeed, here the averaged value (x) takes lower values than xq 
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FIG. 5: Evolution of (x) and ((5x) 2 ) 




related to /? = and the dispersion ((Sx) 2 ) does not decrease in time at late stages. Therefore, the quantity ((Sx) 2 ) 
can be used an effective order parameter. Indeed, if it decreases in time toward zero, then the system moves into 
homogeneous state, whereas nonzeroth values for the dispersion mean that modulated spatial patterns are realized. 
From Figj5] one can see that at f3 ^ at large adsorption rate a the system is in the high density state, whereas 
at elevated interaction strength of the adsorbate small islands of the dense phase are possible; here due to large 
interactions between adsorbate the evaporation/dissolution processes are less probable. Comparing curves related to 
two cases r = and r / one can find that in the last case both (x) and ((Sx) 2 } exhibit nonmonotonic behavior, here 
oscillations with small amplitude are possible. It is well related to results of the linear stability analysis. Moreover, 
it is interesting to note that in the case r = most of islands of the dense phase are of equiaxial symmetry, whereas 
at r 7^ such islands are elongated in one of two possible equivalent directions. 

Let us consider dynamics of the structure function at different values for r and different time intervals (see Figj6|. 
To calculate S(k,t) we have used fast Fourier transformation procedure. Let us start with the simplest case of r = 
(see solid lines in Figsj6^i,b). Here only major peak of S(k,t) is realized that is related to period of islands. There is 
smooth behavior of the structure function tails. During the system evolution the peak is shifted toward stationary 
value of the island size; its height increases (islands become well defined and boundaries between dense and diluted 
phases become less diffusive). In the case r/Owe get one major peak at small wave number and additional peaks 
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FIG. 7: Dynamics of the number of islands (a), averaged island size (b) and the corresponding averaged aspect ratio (c) at 
P — 0.1. Values for the dynamical exponent z related to the scaling law (R) oc t z are shown near the approximate lines. In plot 
(c) snapshot of pieces of the system obtained at t = 200, 220, 240, 280, 300 illustrate a change of R x /R y ; other parameters are: 
e = 4, /? = 0.1, a = 0.1. 



at large k. Emergence of such minor peaks means formation of other patterns (patterns with other periods). In the 
course of time an amplitude of such satellite peaks decreases that means pattern selection processes when the system 
selects one most unstable mode characterized by the major peak whose height increases. This oscillatory behavior of 
the structure function at large wave numbers is well predicted by the linear stability analysis. 

Next, let us study a behavior of the averaged radius of islands. In our computations we have calculated number of 
sites related to one island. This number corresponds to a square of the island. Assuming that an island with spherical 
symmetry has the same square we have computed the radius (R) of the corresponding spherical island. In Fig(7^ we 
plot dynamics of the number of islands at different values for r, e and a. Dynamics of the quantity (R) (measured 
in units of diffusion length Ld) is shown in Figj7]3. It follows that the system attains the stationary state with finite 
number of islands. Comparing curves related to different r, one can find that nonequilibrium effects related to r ^ 
decrease essentially the number of islands (but stationary islands have large size). If the adsoprtion rate a is small, 
then number of islands is large (they are characterized by small sizes). From Figwb it is seen that in the simplest case 
of pure dissipative system characterized by r = the averaged radius is a monotonically increasing quantity. At large 
time interval it attains a stationary value and does not changed. It means that the coarsening procedure is finished 
and we get stationary patterns. It follows that islands of the adsorbate have the averaged length (R) ~ (0.15 -1-0. 9)1/^. 
Following Refs.[34j [35] one can estimate (R) considering deposition of Al on TiN(100): at room temperature one has 
the lattice constant clai = 4.05 x 10 _10 m, the pair interaction energy eo = —0.22 eV with the coordination number 
Z = 4 gives ro = Zclai — 1.6 x 10 _9 m, the diffusion constant D = 10 -10 cm 2 s -1 . Hence the patterns have the size 
(R) ~ (10 -T- 60) x 10- 9 m. 
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x R/<R> ,v R/<R> 

c) d) 

FIG. 8: Evolution of the probability density function of the dimensionless radius of the islands size at different r at f3 — 0.1: 
a) T = 0, e = 4, a = 0.06; b) r = 0, e = 3, a = 0.25; c) r = 0.5, e = 4, a = 0.06; d) r = 0.5, e = 3, a = 0.25 



Studying scaling properties of the island growth we have estimated a scaling exponent z in the islands size growth 
law (R) ex t z . Fitting the data related to stage of growth, we have obtained z = 0.08 for pure dissipative system 
with e = 4. Taking r ^ 0, one can find that at small e (in our case e = 4) the averaged radius of islands increases 
in nonmonotonic manner. Fitting the scaling regime we have found that at r = 0.5 the scaling exponent takes large 
values, in our case z — 0.34. When the adsorption rate decreases the growth processes delay and characterize by 
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small dynamical exponent z (at a = 0.06 one has z = 0.19). Therefore, nonequilibrium effects related to relaxation 
of the diffusion flux characterized by r ^ accelerate growth processes, the same is observed when adsorption rate 
increases. It is interesting to note that at r ^ oscillations predicted for the averaged field and the structure function 
are possible for the dependence (R(t)}. It means that islands can change their size during the system evolution in 
oscillating manner. Amplitude of oscillations is well pronounced for large a and small e. 

To prove that islands can oscillatory change their sizes we compute an averaged aspect ratio (AR) = (R x /R y ), 
where R x and R y are sizes of an island in x- and ^-direction, respectively. From the obtained dependencies shown in 
FigjTJ^ it follows that for pure dissipative system deviations from the straight line are small and may be considered as 
fluctuations in R x and R x values. At r ^ and small e oscillations in (AR) are well pronounced. Such oscillations 
in lateral and longitudinal sizes of islands mean that at some fixed time interval most of the islands grow in one 
direction whereas in other direction their size decreases, in next time interval these two directions are changed. When 
this scenario is repeated one gets an oscillation picture of island size growth. In FigjTf; we present snapshots showing 
change of the growth orientation during the system evolution. 

In our study we are interesting also in a behavior of the probability density functions (PDF) of the island sizes 
distribution. Let us consider PDF's for pure dissipative system, initially. As Figs(8^i,b shows islands of different sizes 
formed at initial stages evolves in such manner that the most probable value for R is located around (R). Islands of 
other sizes are possible with small probabilities. In the stationary limit PDF has unique well pronounced peak, small 
fluctuations of island size distribution are possible and are shown as satellite peaks around major one. In the case 
of r 7^ (see Figsj^d) the hyperbolic transport promotes formation of islands with different sizes. The principle 
difference comparing to the case of pure dissipative system is in a bimodal form of PDF observed at large times. Here 
at early stages islands of different sizes are formed, and during the system evolution following the Ostwald ripening 
mechanism small islands dissolve and large islands reduces their sizes. In such a case two well pronounced major 
peaks of PDF indicate that there two most probable sizes of islands related to R > (R) and R < (R) . 

V. CONCLUSIONS 

We have studied dynamics of islands formation of the adsorbate using generalized approach including persistent 
motion of particles having finite speed at initial stages and diffusion kinetics at final ones. It was found that stabi- 
lization of nano-patterns in such class of reaction- Cattaneo models is achieved by noneqilibrium chemical reactions. 
It was found that during the system evolution pattern selection processes are realized. We have shown that possible 
oscillatory regimes for islands formation are realized at finite propagation speed related to nonzero relaxation time 
for the diffusion flux. 

Our results can be used to describe formation of nano-islands at processes of condensation from the gaseous phase. 
Despite we have considered a general model where relaxation time tj for the diffusion flux is small but nonzeroth value, 
one can say that condensation processes with formation of metallic islands can be described in the limit tj/lj^ 1 < 10 -3 
(here ujd is the Debye frequency), whereas nano-islands formation with tj/uj^ 1 ~ 10 _1 10 -2 is possible for soft 
matter condensation (semiconductors, polymers, etc.). 
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